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Abstract 

We propose two main applications of Gyongy (1986)'s construction of 
inhomogeneous Markovian stochastic differential equations that miniick 
the one- dimensional marginals of continuous Ito processes. Firstly, we 
prove Dupire (1994) and Derman and Kani (1994)'s result. We then 
present Bessel-based stochastic volatility models in which this relation is 
used to compute analytical formulas for the local volatility. Secondly, we 
use these mimicking techniques to extend the well-known local volatility 
results to a stochastic interest rates framework. 



t:::j- ■ 1 Introduction 

O' 

^O ■ It has been widely accepted for at least a decade that the option pricing theory of 

^p I Black and Scholes (1973) and Merton (1973) has been inconsistent with option 

r^ ■ prices. Actually, the model implies that the informational content of the option 

surface is one dimensional which means that one could construct the prices of 
r^ ■ options at all strikes and maturities from the price of any single option. It has 

J~J I also been shown that unconditional returns show excess kurtosis and skewness 

^ ■ which are inconsistent with normality. Special attention was given to implied 

K^ [ volatility smile or skew, but research has concentrated on implied Black and 

j_j ■ Scholes volatility since it has become the unique way to price vanilla options. 

C^ I Accordingly, option prices are often quoted by their implied volatility. Nev- 

ertheless, this method is unsuitable for more complicated exotic options and 
options with early exercise features. To explain in a self-consistent way why 
options with different strikes and maturities have different implied volatilities 
or what one calls the volatility smile, one could use stochastic volatility models 
(eg. Heston (1993) or HuU and White (1987)) 

Given the computational complexity of stochastic volatility models and the 
difficulty of fitting their parameters to the market prices of vanilla options, 
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practitioners found a simpler way to price exotic options consistently with the 
volatility smile by using local volatility models as introduced by Dupire (1994) 
and Derman and Kani (1994). Local volatility models have the advantage to 
fit the implied volatility surface; hence, when pricing an exotic option, one feels 
comfortable hedging through the stock and vanilla options markets. 

In the last twenty years, academics and practitioners have been primarily 
interested in building models that describe well the behavior of an asset whether 
it is equity, FX, Credit, Fixed-Income or Commodities and very rarely models 
that specify any cross-asset dependency. For all cross-asset derivative products, 
this dependency modifies the model one should use or at least the calibration 
procedure. Certainly, models that incorporate a dependency on other asset 
classes than a specific underlying need to be recalibrated as soon as the other 
asset classes become random, in particular in the fast growing hybrid industry 
where it is necessary to model several assets. 

The remainder of the paper is organized as follows. Section 2 recalls prelim- 
inary results on Bessel processes and states mimicking properties of continuous 
Ito processes exhibited by Gyongy (1986) and Krylov (1985). Section 3 re- 
calls well-known results of Dupire (1994) and Derman and Kani (1994) on local 
volatility, gives a proof of the existence of a local volatility model that mimicks 
a stochastic volatility one based on Gyongy (1986) theorem. Section 4 gives ex- 
amples of stochastic volatility models where a local volatility can be computed. 
Those examples are based on remarkable properties of Bessel processes such 
as scaling properties. In order to extend the class of volatility models (where 
closed-form formulas can be obtained), we propose a general framework in which 
the volatility diffusion is a general deterministic time and space transformation 
of Bessel processes. Analytical computations are proposed in cases where the 
volatility diffusion is independent from the stock price diffusion as well as in 
cases where they are correlated. Section 5 applies the results of Section 3 to the 
case of stochastic interest rates and more generally shows how Gyongy (1986) 
theorem can be applied to construct a local volatility model in a deterministic 
interest rate framework, starting from a stochastic volatility model with stochas- 
tic rates. Finally, Section 6 concludes our work and presents an important open 
question on mimicking the laws of Ito processes. 

2 Preliminary Mathematical Results 

2.1 Bessel and CIR Processes 

Let {Rt,t > 0) denote a Bessel process with dimension 6, starting from and 
iPt,t > 0) an independent brownian motion from {Bt,t > 0) the driving brow- 
nian motion. Let us recall that Rj solves the following SDE: 

di?^ = 2RtdBt + Sdt 



and let us now define 



RsdPs 



and At 



Rids 



Then, the one-dimensional marginals of (At, It) are at least in theory well- 
identified, via Fourier-Laplace expressions, and are closely related with the so- 
called Levy area formula (see Levy (1950), Wilhams (1976), Gaveau (1977), Yor 
(1980), Chapter 2 of Yor (1992) and many other references). Here we simply 
recall, for our purposes the formulae: 
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as well as: 
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a formula that we shall use later. Some developments for the law of At are 

given, e.g. in Pitman and Yor (2003). 

For a Bessel process of dimension 8 starting at x, one gets the following formula: 
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Let us now present a scaling property of the Bessel process with respect to 
conditioning, which is important in the sequel. 



Proposition 2.1 For any Bessel process Rt with dimension 5, we have: 
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(3) 



Remark 2.2 This result is in fact a very particular case of a more general 
result involving only the scaling property of the process {Rt,t > 0), see, e.g, 
Pitman and Yor (2003). But, for the sake of completeness, we shall give a 
direct proof of \^ below: 

Proof : From the scaling property of {Rt,t > 0), we deduce that for every 
/ G C^(R, R+), with bounded derivative, we have: 
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We then differentiate both sides with respect to i, to obtain: 
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Since this identity is true for every bounded Borel function /', the identity ((3J) 
foUows. ■ 

Remark 2.3 We now check that formula ^ can be obtained directly as a con- 
sequence of formula 0); differentiating 0) both sides with respect to a and 
taking a = 0, we obtain: 
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while, taking a ^ in \^, and differentiating both sides with respect to b, we 
get: 



bE 



Atexp{-jAt) 



St 



Mbt) 



2(cosh(6i))2+i 
and the identity 0) follows from the comparison of these last two equations. 

A reason why squared Bessel processes play an important role in financial 
mathematics is that they are connected to models used in finance. One of these 
models is the Cox, Ingersoll and Ross (1985) CIR family of diffusions which are 
solutions of the following kind of SDEs: 



dXt 



bXt)dt + a^\x7\dWt 



(4) 



with Xq — xq > 0, a E IR+, b E M., a > and Wt a standard brownian motion. 
This equation admits a unique strong (that is to say adapted to the natural 
filtration of Wt) solution that takes values in M+. 

One is now interested in the representation of a CIR process in terms of a 
time-space transformation of a Bessel process: 

Lemma 2.4 A CIR Process Xt solution of equation Op can be represented in 
the following form: 



Xt 
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where R denotes a Bessel process starting from xq at time t = of dimension 
- -^ 

Proof : This lemma results from simple properties of squared Bessel processes 
that can be fomid in Revuz and Yor (2001), Pitman and Yor (1980, 1982). ■ 

This relation has been widely used in finance, for instance in Geman and 
Yor (1993) or Delbaen and Shirakawa (1996). 



2.2 Mimicking Theorems 

A common topic of interest of Krylov and Gyongy respectively in Krylov (1985) 
and Gyongy (1986) is the construction of stochastic differential equations whose 
solutions mimick certain features of the solutions of Ito processes. The con- 
struction of Markov martingales that have specified marginals was studied by 
Madan and Yor (2002). Bibby, Skovgaard and S0rensen (2005) as well as Bibby 
and S0rensen (1995) proposed construction of diffusion-type models with given 
marginals. 

Let us now consider an Ito differential equation of the form: 
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where Wt is a JFi-Brownian motion of dimension k, ((5t)tgR^ and (/3t)tgR^ are 
bounded JF^-adapted processes that belong respectively to M„jt(R), the space 
of n X fc real matrices and to M" . 

Definition 2.5 (Green Measure) Considering two stochastic processes Xt, 
valued in M" and jt, with jt > 0, one defines the Green measure fJ-x.-y by: 



Aix,7(r) 
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dt 
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where T is any borel set ofW^ 

Remark 2.6 The stochastic process 74 is called the killing rate 

Theorem 2.7 (Krylov) If ^t is an ltd process defined as previously and sat- 
isfying the uniform ellipticity condition: 3X G K+ such as SS* > XIn 
as well as the lower houndedness condition: 3a G M+ such as j > a, 
then there exist deterministic functions a : M" — > 
c : R" ^ R+ such that the following SDE: 



dxt 
has a weak solution satisfying: 

vr e B{i 

fOO 

E 



= <T{xt)dWt 

= 



" ^ M„(M) 
b{xt)dt 



lr(a;t)exp(— / c{xs)ds)dt 



lr(6)exp( - / jsds)dt 
Jo 

ie : M«,7(r) = A*x,c(r) 

Proof : See Krylov (1985) ■ 

Definition 2.8 (Weak Solution) The stochastic differential equation 

dXt = f{t,Xt)dWt+git,Xt)dt (8) 

Xo = (9) 



is said to have a weak solution if there exist a probability space (fi, J^, P) and an 
Tt-Brownian motion with respect to which there exists a J^t- adapted stochastic 
process Xt that satisfies 0) and 1^). 

A natural question asked and answered by Gyongy is whether it is possi- 
ble to find the solution of an SDE with the same one-dimensional marginal 
distributions as an Ito process. The answer is stated below: 

Theorem 2.9 (Gyongy) If ^t is an Ito process satisfying the uniform elliptic- 
ity condition: 3X e W^ such as 66* > XIn 

then there exist bounded measurable functions a : M+ x M" -^ M„^„(M) and 
6 : R+ X M" ^ R defined by: 

V(t, x) eR+x M" 









a{t,x) 


= ^6t6;\^t^: 

V 








b{t,x) 


- ¥.[Pt\it = x 


such that the 


following 


SDE: 










dxt 


- C7(i, 


Xt)dWt +b{t,xt)dt 
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has a weak solution with the same one- dimensional marginals as ^. 
Proof : See Gyongy (1986) ■ 

Remark 2.10 Two kinds of mimicking features of a general ltd process were 
illustrated in this section . With Krylov, we were able to construct a Markov 
homogeneous process solution of an SDE, that has the same Green measure 
than the ltd process. Using Gyongy 's results, we were able to build a time- 
inhomogeneous Markov process solution of an SDE that has the same one- 
dimensional marginals as the Ito process. 

A possible extension to the above mimicking property is to consider a real 
Ito process ^ driven by a multidimensional Brownian motion and obtain a new 
mimicking result useful for the remainder of the paper; the proof is straightfor- 
ward from Gyongy (1986) proof. Let ^ be as follows : 

6= / <6„dW,>+ f Ms (10) 

Ja Jo 

where Wt is a .7-"t-Brownian motion of dimension k, {6t)tem.^ and (/3t)tgR^ are 
bounded JFt-adapted processes that belong respectively to M'^ and to M. 



Theorem 2.11 If ^t 'is cin ltd process defined as in i^lUJI satisfying the uniform, 
ellipticity condition: 3X g ]R!j_ such as 66* > A 

then there exist bounded measurable functions a : M+ x R ^ M and b : M+ x M — > 
R defined by: 



V(i,x) eR+ xR 




a{t,x) = 


-- (E[6t6;\^t-^ 


b{t,x) = 


= E[pt\^t=x 



such that the following SDE: 

dxt = o-{t,xt)dWt +b{t,xt)dt 

2^0 = 
has a weak solution with the same one- dimensional marginals as ^. 

3 Generalities on Local Volatility 

3.1 Fokker-Planck Equation 

Let us assume that: 

^^r{t)dt + a{t,St)dWt (11) 

where r and a are deterministic functions, a is usually called the local volatility. 
Under the local volatility dynamics, option prices satisfy the following PDE: 

and terminal condition V{S,T) = C{S,T) = PayOffriS). 
If we consider call options, we would get V{S,T) — {S — K)+. It has been 
proved that one can obtain a forward PDE for C{K, T) instead of fixing {K, T) 
and obtaining a backward PDE for C{S, t). To get the Forward PDE equation, 
one could just differentiate H12|l twice with respect to the strike K and then get 
the same PDE, with variable (f> — ^^ and terminal condition 6{S — K). (p is 
the transition density of S and is also the Green function of (|12|l . It follows that 
as a function of {K,T) satisfies the Fokker-Planck PDE: 

Now, integrate twice this equation taking into account the boundary conditions, 
one obtains the Forward Parabolic PDE equation: 

dC aHT,K) d^C dC_ 

dT 2 ^ aif2+nJ^A^^-u {16) 



with initial condition C{K,0) = (Sq — K)-^-. Hence, one obtains Dupire (1994) 
equation 

-^m^) = ^4^iW^ (14) 



Moreover, if one expresses the option price as a function of the forward price, 
one would write a simpler expression: 

dc 
aHT,K,So]- '^ 



2" dK^ 



where C is now a function of (-Ft, K, T) with Ft = Sq exp ( J„ r{s)dsj . 

3.2 Matching Local and Stochastic Volatihties 

A stock price diffusion with a stochastic volatility is one of the following form 

dSt 



^ = r{t)dt + ^VtdWt (15) 

where Vt is a stochastic process, solution of an SDE and r{t) is a deterministic 
function of time. (We do not yet discuss the dependence of the stock price and 
volatility processes, also called Leverage effect) 

One can find a relation between the local volatility and a stochastic volatility. 
First, one applies Tanaka's formula to the stock price process: 

e-/oK^)'^-(5,_x)+ = {Sq-x)+- r{u)e-^«^'^'^'''{Su-x)+du 




t -1 i-t 
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Assuming that (e '0 '"^"'■^"S'i, i > 0) is a true martingale, then 
(/„ l[Su>x}d{e^ '0 ^^^^'^^Su),t > 0) is a martingale and one gets: 

E[e--fo'^(^)'^^(5,-a;)+] = E[{So - x)+] + x f E[r{u)e- fo '<''>^n{s^>,}]du 

Jo 

If/"* 1 

+ E / e-'C '-('')'''' di^ (5) 

Then, differentiating the previous relation and using Fubini theorem, one ob- 
tains: 

dtC{t,x) = xE[r{t)e--f'><'^^n^s,>.}]dt+^ne-'^''^^''>'''dmS)] (16) 

where C{t, x) — E[e^'o '"("''^''(54 — a;)+] Using a classical characterization of the 
local time of any continuous semi-martingale: 



1 /•* 
mS) = lim - / ^{x<s,<x+e}d < S,S >s 



(17) 



one gets with a permutation of the differentiation and the expectation: 

(18) 
as a result of d < 5,5 >t— VtSfdt. Now, one may write using conditional 
expectations and the fact that interest rates are assumed to be deterministic, 
the following identity: 

E[VtSfl{,<s,^^+,}] = E[E[14|5t]52l{,<s^<,+,}] 

From there, one easily obtains: 

\im-E[VtSfl[^<s,<..+e}] = limiE[E[l^t|5t]5,2l{,<s,<.+.}] 

= E[Vt\St=x]x''qtix) 

where qt (x) is the value of the density of St in x. Since Breeden and Litzenberger 
(1978), it is well known that ^^|^ = e' -fo '■(^)d.^qt{x). It is also known that 
H = -E[e-fo'''-''^'^n{St>x}] One finally may write: 

f + „«)f .EKIS.^4.^0 (19) 

Comparing equation H14|l and the above equation, one may obtain an equation 
that relates local and stochastic volatility models 

a^{t,x)=E[Vt\St=x] (20) 

Hence, we have proven that if there exists a local volatility such as the one- 
dimensional marginals of the stock price with the implied diffusion are the same 
as the ones of the stock price with the stochastic volatility, then the local volatil- 
ity satisfies equation ^^. 

Another way to prove this relation is to apply Gyongy (1986) result. Since 
the stock price dynamics with a stochastic volatility given by equation H15|l and 
the ones with the local volatility given by equation Hll|l must have the same one- 
dimensional marginals, one can apply Gyongy Theorem: assuming that there 
exists A e M![. such that 5^^ > A we get the well-known relation between the 
local and the stochastic volatilities: 

a{t,St^x)^ (E[Vt\St^i 

It is important to notice that Gyongy gives us the existence of such a diffusion 
in addition to provide an explicit way to construct it. More generally, assuming 
just that the volatility process is a general continuous semi-martingale, one can 
also get the same result, and a justification for the use of local volatility models. 
Hence, we obtain an illustration of Gyongy's result in a finance framework. 



Moreover, it is shown that one can get the relation H2U|I without using the 
Forward PDE equation. 

As a first remark, we should notice that if we choose Vt such as ^/Vt = 
(7{t, St), we then obtain another direct proof of equation (|14|l . 

As a second remark, we can prove that if {St — e^ ■>« ''*^*''**S't), t > 0) is a 
strict local martingale (which is studied in Cox and Hobson (2005) who named 
this market situation a bubble), then 

nl Ms^>.}dSu] - nst-s^]~nl Ms^<.}dSu] 

Jo Jo 

= nst-so] 

since using Madan and Yor (2006), (/g l{s„<^}d(e-^" ''('')'*'' S'„),t > 0) is a 
square integrable martingale. Hence, defining 

cg(i) = E[So - St] 

and assuming that cs is a continuously differentiable function, one obtains an 
extension of equation (|14|l that is a generalization to the case of strict local 
martingales. This equation writes 



f +a;r(i)f +c^(t) 



2/, \ _ dt ' "^ ' Ox ' S^ 
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4 Applications to the Heston (1993) model and 
Extensions 

4.1 The Simplest Heston Model 

The aim of this paragraph is now to compute the local volatility not by excerpt- 
ing it from the option prices (see for instance Derman and Kani (1994)) but by 
applying Gyongy's theorem. 

Among the possible choices of stochastic volatility models, we will consider 
the simplest one, given by the following SDE: 

^ = WtdBt (21) 

So^l 

where {Wt) and {Bt) are two independent one-dimensional Brownian motions 
starting at 0. We do not consider any drift term in our stock diffusion as we 
look at the forward price dynamics that are driftless by construction. 

To make our discussion a little more general than the model presented in 
equation l|^ . we write l(5lj: 

^ = \Wt\sgn{Wt)dBt 
10 



Now we define f3t = /g sgn{Ws)dBs, another Brownian motion which is in- 
dependent of {Wt,t > 0) and consequently of the reflecting Brownian motion 
(|l^t|,i > 0). We get the foUowing model: 

7 Q 

Now this modified form leads itself naturally to the generalization: 

^ = Rtd/3u So = 1 (22) 

where, as in subsection 2.1, {Rt) denotes a Bessel process with dimension 5 
starting at and (/3t) an independent Brownian motion. 

Let us consider a Markovian martingale (St,i > 0), which is the unique 
solution of: 

^=a(i,I],m (23) 

for some particular diffusion coefficient {a{t,x),t > 0,x G M+j which has the 
same one-dimensional marginal distributions as {St,t > 0) the solution of (|22(l . 
We will now use DroDOsition l2.1l to find a, the local volatility. Wc follow the 
notation in subsection 2.1, and introduce a useful notation: 



N^At-fiAt (25) 



L\^> = It-fiAt (24) 

(Law) 



where iV is a standard gaussian variable independent of At . Next we remark as 
a consequence of (|25|) that for any fixed i > : 

(i?„Ll^V'^-"'^i?*,^V^-M^) 

and 

E[i?2|L('') =l]= E[EiR''t\N,At)\Ny/At ~ fiAt = /] 

Since N is independent of Rt , we obtain 

E[i?2|L(^' =;] =E[E{R^t\At)\N^t-f^At^l] 
From Q, we deduce: 

E[i?2|L(^) = ?] ={^)E[At\N^t~f^At = l] (26) 

Now, the computation of the expression in H26|) is a simple exercise, which we 
present in the following form: 
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Lemma 4.1 Let X > be a random variable independent from a standard 

gaussian variable N . Denote Y^^' = NyX — fiX . Then: 

i) for any f : M+ — > M+ , Borel function, the following formula holds: 



E[/(x)|r(^) = z]- 



/l(*')(l;z) 






2X 



where: h^^'{f; z) - 

ii) in particular, for f{x) — x, one can write . 

/ dk 
= -I Ik 

k 



:[x|y(^) 



z2 ,,2 



(27) 



where k{a, 6) = E 



1 exp -(f +6X) 



The proof of this lemma resuhs from elementary properties of conditioning 
and is left to the reader. 

We now give a formula for a'^(t,x) in terms of the law oi At = A]. ' , by using 
equation lPB)l and the above lemma. Indeed, it follows from these results that: 

^ ab V 2 ' 8/ 



E[R^\\niSt) = I] 



* ^Kt's) 



where fc^(a, 5) = E -^ exp - (^ + bAt) 

Using the scaling property, we have kg{a,b) = jfc](-p-,6i^) which allows us 
to concentrate on k^{a, b) = kg{a, b). 

The following formula for the density fs of Ai is borrowed from Biane, Pitman 
and Yor (2001). Denoting h — ^, we have: 



fsix) EE flix) 



T{h)2-^ ') r(n + l) V2^ '"P 



{2n + hy 
2x 



(28) 



For 6 = 2, A^'^' , or equivalently f2{x) = /{(a;) enjoys a symmetry property (also 
shown in Biane, Pitman and Yor (2001)): 
For any non-negative measurable function g 
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(29) 
(30) 

(31) 



n=0 



12 



From formula (|28|l . one may compute with the change of variables a — ^,b 



2 /o2 



k\a,b) = E 



' -pf-^(^+/3^^'^0 



Vlw 



2M(^ 
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7I)£^j" ^ r(n+l) V2^ io ^' ' 



-+/3^2:) 
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Also of importance for us, is the result: 
(h^(n hW = — ^ 

dh 






oOg-i(a2x+^ 



r(/.) 



n=0 



X 



where a„ = ^a^ + (2n + /i)^. 

Recall the integral representation for the Mc Donald functions K^: 



1 ,z. 



dt 



1 ^^\-{.^^4) 



In particular, we have: 

Ko{z) 
As a consequence: 

r^e-^("^-^)=2i^of^^ 
Now, plugging (|34|l in (|32|l . we obtain: 



2./o i 



^fc^(a,6))- ^^'^ 
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r(/^) 



Likewise, we deduce from (|55|) that: 



which implies 



Ki{z) = K_i{z) ^ - dte-y^^"^' 
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Hence, we get as a consequence of (• 

2h °o 
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T{h) ^ 

^ ' n=0 
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(33) 



(34) 



^ r(n + /.) 2n + /. /g^X 

^^ ^ rn + 1 72^ %V2y ^ ' 

n=0 \ y V V 



(36) 
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Recalling that (3 — v26 and that q;„ = \/2a + (2n + h)'^, we may now write 
^ and ^ as: 

X, , 2'=- ^, ,„Tin + h)2n + hV2b , rr, 

^ ^ n— ^ 

And finally, we obtain the following formula for the local volatility: 

A / W^\ /2 +2 
"^('•--''-l(^)<i-^' (39) 

4.2 Adding the Correlation 

We now assume a non-zero correlation between the volatility process and the 
stock price process. This is a common fact in finance called the Leverage Effect 
and translated by a negative correlation. For a financial understanding of this 
effect, one can refer for instance to Black (1976), Christie (1982) or Schwert 
(1989). 

Let us define our new model for the stock price dynamics with a Bessel 
process of dimension 5 starting from correlated to the Brownian motion of the 
stock price process: 

^ = RtdWt (40) 

dRJ = 2RtdWt + 5t (41) 

d<W",W >t= pdt (42) 

Sq = 1 and Ro = (43) 

Then, there exists a Brownian motion (3 independent of the Bessel process 

such that Vt: 

Wt = pW^ + Vi - p'Pt 

Using the previous formula, plugging it in (|40|l and then inserting H41|) in the 
new (pn| . one gets: 

^ = ^ {dRl - 5dt) + ./l~^RtdPt (44) 

Then using Ito formula applied to /(x) — ln(x) 

dSt 1 ^2. 



d\n{St) ^ ^ - -R^dt 



one obtains: 



In(St) = f (^? - St) + V'^-P^J Rsdps ~\j Rids 



(45) 
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Let us consider as in subsection 5.1, L^' — J^ Rgdfis — A^/q R^s'^s (we are 

= ). Since R and [3 are independent, 



especially interested in the case /^ = — j^ 

we shall use the same notation as above. Particularly, At and It will refer to 

the quantities defined in subsection 2.1. 

Now, equation (|45|l can be rewritten as follows: 



HSt) = f (i?? - 5t) + ^l^Lt'^^^' 



(46) 



Since we wish to evaluate the local volatility E[i?f | ln(5f) = Z] , we will try to 
compute more generally the following quantity: 



¥.[Rl\mRl + U^^' ^i\ 



(47) 



where m is a real constant. 



Remark 4.2 We immediately see that if we take to = 0, i.e p = 0, we are back 
to the previous paragraph setting. 

First, we see that equation 1)26(1 is easily extended to the case with correlation 
and we obtain: 



(48) 



E[R'^\mR'^ + L\^> =1] = -ElAtlmR^ + L\^> = l] 
Before extending Lemma 4.1, one must recall that for any i > 0: 

iRt,At, l[^^) ^''^"'^ (i?„ A, N^t ~ M) 

where iV is a standard gaussian variable independent of Rt and At . The follow- 
ing simple result will be helpful for the remaining of the paper: 



Lemma 4.3 Let X > and Z > independent from a standard gaussian 

variable N. Denote F^^) == Ny/X-pX. Then: 

i) for any Borel function f : R^ — > R+ , real number m we have the formula: 



E[f{X,Z)\mZ + Y^^'^ 



-(ti.,m) 



(/;^) 



a(^'™)(l;2) 



where: a''^'•"'\f;z) =E 



f{X,Z) 



exp 



2X 



a) in particular, for f{x,y) = x, we obtain 



M\/2l a MV2'V2'V2^ 



where a{a, b,c) —E, 



Vx 



exp 



X 



9x 



bcZ)] 



(49) 
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The other fundamental resuh we now need, is the joint density of {R^, L R'^ds)t>o. 
Theorem 4.4 The joint distribution gt of {R^,L R^ds) is given by: 



g,{x,y) - ^ , X]^-^'^^'V^'^'Vi'(^,y) (50) 



where fl is defined by 

Dv{C) is a parabolic cylinder function and {v)k the Pochhammer's symbol defined 
by iiy)k = y{v + l)...(z/ + fc - 1) = r(i. + k)/V{v) 

Proof : See Ghomrasni (2004) who evaluates the Laplace transform of 101 in 
order to get the density function. ■ 

For the definition and properties on the parabolic cylinder functions, we refer 
to Gradshteyn and Ryzhik (2000). 
Let us define at in the following form: 



at{a^ 6, c) = E 



1 ( Aa- cR^^^ 



exp ( - (^ P^ + b'^At + bcRl" 



(52) 



.V^t "^K ' At 

Unfortunately, there is no more scaling property as in the zero-correlation case 
and we may not rewrite a* as a function of t and ai . One can then compute 
the local volatility a^''' by noticing that in the case of particular interest for us, 
the parameters are defined as follows: 

"^ — and z = — and ji = 



2v/r^7 vi^p^ 2v/r~o^ 



We then obtain 

a^p\t,x^e')^-^2{l-p^) 



'at \ 1 M EXi 



at y V2(l-P^)' V8(l-P')' VW^T)' 

(53) 



4.3 Prom a Bessel Volatility process to the Heston Model 

The Heston (1993) model for representing a stochastic volatility process is a 
particular case of the Cox, IngersoU and Ross (1985) stochastic process, of the 
form: 

dVt = k{9 - Vt)dt + rjy/vidWt (54) 

with initial condition Vq — vq 

Actually, it is possible to find out deterministic space and time changes such 
as the law of the Heston SDE solution and the Time-Space transformed Bessel 
Process are the same. 
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Proposition 4.5 For every Heston SDE solution, there exist a Bessel process 
and two deterministic functions f and g with g increasing such as: 

Vt = fit) X i?2^,) 

where R denotes a Bessel Process of dimension S — ^^ starting from ^/vq at 
time t = and f and g are defined by: 

fit) = e-«* 

git)^^ie-'^l) 
Proof : It is just an application of Lemma 2.4. ■ 

One may now apply the results of the previous sections using the time and 
space transformations presented in the previous paragraph 

Proposition 4.6 Let us consider the following stochastic volatility model: 

-— =; ^/v^dPt, S{t=o} == So 

dVt = K{e - Vt)dt + rjy^tdWt 
d< P.,W. >t = pdt 

where Pt is a Brownian motion and Vt is an Heston process as defined above. 
Then the local volatility a that gives us the expected mimicking properties, sat- 
isfies the following equation: 

a(t,.) = ^e-a(|^(e--l),|-) (55) 

where a'^{t,x) = IE[i?^| exp(/f — hAt) ~ x] and Rt is a Bessel Process of dimen- 
sion 5 — -^ starting from Vq . 

Proof : First, one has the Gyongy volatility formula: 

a^t,x)^^e^-'E[Vt\St=x] (56) 

Then using Lemma 2.4, one easily obtains the result. ■ 

Remark 4.7 Let us note that we only have closed-form formulas in cases where 
Vo = and that otherwise we have to go through Laplace transform inversion 
techniques. 

One can propose a general framework for constructing stochastic volatility 
models based on Bessel processes. Local volatilities can be computed through 
the proposition below whose proof is left to the reader. 
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Proposition 4.8 Let us consider the following stochastic volatility model: 
dSt 






dVt - 
d< P.,W. >f = pdt 



5f{t)9'{t) + J^Vt)dt + y^g\t)^tdWt 



where j3t and Wt are Brownian motions, f is a positive continuously differen- 
tiable function and g an increasing C^ function. 

Then the local volatility a that gives us the expected mimicking properties, sat- 
isfies the following equation: 

a(i,x)=5'(t)c7(.g(t),^) (57) 

where a'^{t,x) = E[i?f | exp(/( — ^At) — x] and Rt is a Bessel Process of dimen- 
sion S starting from Vq . 

5 Pricing Equity Derivatives under Stochastic 
Interest Rates 

5.1 A Local Volatility Framework 

With the growth of hybrid products, it has been necessary to take properly 
into account the stochasticity of interest rates in FX or Equity models in a way 
that makes the equity volatility surface calibration easy at a given interest rate 
parametrization. It has been now a while that people have been considering 
interest rates as stochastic for long-dated Equity or FX options, but they have 
not been thinking about it in terms of calibration issues. Besides, according to 
the interest rates part of an equity - interest rates hybrid product for example, 
the instruments on which the interest rates model will be calibrated are different; 
hence it becomes necessary to parameterize the volatility surface efficiently. For 
most of hybrid products, no forward volatility dependence is involved and then 
a local volatility framework is sufficient. Let us now consider a local volatility 
model with stochastic interest rates: 

^ = rtdt + ait,St)dWt 

where rt is a stochastic process and a a deterministic function. 

Now, we can observe that equation H18|l is still valid under stochastic rates and 

we may then write 

dtCit,x) = xE[rte-/o '■'■'^n{s^>,}]rfi+i lim E[il{,<s,<,+,}e-^''-»''V2(i, SJ^^j^^ 
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The second term of the right-hand side may be written as foUows 
and then we have 

where qt{x) is the value of the density of St in x. It is easily shown as well that 



^qt{x)E[e-fo-^ds^g^^^^ 
ux- 

' define the t-fon 
(1989)) by 



dx'^ 
Let us now define the t-forward measure Q* (see Geman (1989), Jamshidian 



-^ = —, where B(0, t) = Efe" ■''o '■=''*1 

dQ B{0,t) V , ^ L J 

Hence, we finally obtain an extension of Dupire (1994) 's formula : 



2(, ,_ ^-xB{0,t)e[ni[s.^,}] 

2 'dx^ 

Under a T-forward measure for T >t, one has 

IE^[rT|^t]-/(t,T) 

where f{t,T) is the instantaneous forward rate. To conclude this subsection, 
we can first notice that this slight extension of Dupire equation may be also 
written 

2 l)x^ 

We then assume that it is possible to extract from markets prices the quantities 
Cow*(rt; lr-g^^j.1 ) (i.e. there exist tradeable assets from which we could obtain 
these covariances) in order to add stochastic interest rates to the usual local 
volatility framework. For the remainder of the paper, we denote this assumption 
the (iJC)-Hypothesis that stands for Hybrid Correlation hypothesis. Under 
this market hypothesis, one is able to calibrate a local volatility surface with 
stochastic interest rates implied by the derivatives' market prices. 

5.2 Mimicking Stochastic Volatility Models 

In this subsection, we consider the case of a stochastic volatility model with 
stochastic interest rates and see how it is possible to connect it to a local volatil- 
ity framework. Let us consider the following diffusion 

dSt 



ndt + ^VtdWt 
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with Vt a stochastic process and let us use equation (I18f) in order to exhibit a 
new mimicking property: 

Then, 



x^ 



d^CE[Vte-Io'-^'^'\St^x] 
dx^ E[e-fo'-^'^''\St=x] 



Hence, we obtain 

E[Vte-fo-^ds^St=x] _ ^ + xf{0,t)^ - xB{0,t)Cov\rt;l{s,>.}) 



]E[g-/„'r,ds|5.^^ 



2 'dx^ 



Spot Mimicking Property Finally, if there exists a stochastic process, so- 
lution of the following SDE 

^ = rtdt + a{t,Xt)dWt 

such that the one-dimensional marginals of the triple ('"ti/g i"sds,Xt) are the 
same as (rj, L rsds, St), then by identification one must have 

2,^ , E[yte-/o'-=''^|5t=x] 



a^{t,x) 



E[e-fo'-^'^'\St^x] 



The existence is easily proven in the cases where {rt,t > 0) is a Markovian 
diffusion. Hence, we exhibit a strong mimicking property since we obtained an 
explicit way to construct a local volatility surface. 

Remark 5.1 We may notice that if interest rates are deterministic, we recover 
the well-known formula \2U\) . 

Forward Mimicking Property Let us now write a Forward mimicking prop- 
erty by applying Gyongy's result to match the one dimensional marginals of a 
stochastic volatility model and of a local volatility one: 

If one defines F^^^ = Stc' ^o"^''^" and F^"^^ = Xte' ^o^''^" where S and X are 
defined above, we obtain the existence of diffusions Y"/ and Y^ solutions of 



dY, 



(«) 



* — V, . /'+ v(i) 



yi^) 



= ^,)itX')dWt 
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for i = 1,2 such as 

Since the onc-dimcnsional marginals of Fj and Ff; must be equal, one obtains 

E[Vt\St = xef» ''^'^'] = E[cr2(t, xe- -fo ''^'^')\Xt = xe-^o ''=''''] (59) 

We consequently obtain an implicit way to construct a local volatility surface we 
say that this relation is weak in the sense that it is a weak mimicking distribution 
property which is involved in the above relation. 

5.3 Prom a Deterministic Interest Rates Framework to a 
Stochastic one 

Going from a framework to another is valuable for calibration issues. Let us 
assume, for instance that a model has been calibrated with deterministic inter- 
est rates and that one wants to recalibrate the same model assuming stochastic 
interest rates. Let us introduce some notation to define the different kinds of 
frameworks we will go through in this subsection. 

Notation 

LV stands for Local Volatility, SV stands for Stochastic Volatility, DIR stands 
for Deterministic Interest Rates and SIR stands for Stochastic Interest Rates 

From DIR-LV to SIR-LV Let us first consider the local volatility case. 
Under deterministic interest rates, the stock price dynamics are driven by the 
equation 

^^r{t)dt + a{t,St)dWt 
while under stochastic interest rates it would be 

^ = ndt + a{t,St)dWt 

and we know that both local volatility functions solve the following implied 
equations: 



a^{t,x) 
a'^{t,x) 
where /(0,t) = r(f). 



2 'dx^ 



^ + xf{0,t)^ - xB{0,t)Cov\rt;l{s,>.}) 



2 'dx^ 
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Now, if the prices involved in the estimation of the local volatility surfaces 
are observed on markets and respect the (_ffC)-Hypothesis, one may write 

a {t, x)-<7 {t, x) ^ -^^ (60) 



Prom DIR-SV to SIR-SV If we assume that a general Ito process drives 
the volatility we will write 





= T(t)dt + ^V^^^dWt 


dSf^ 


= rtdt + \lvi:'^UWt 



and then, if S^^' and S^"^' have the same one-dimensional marginals, we obtain 
the following relation to relate V'^^^ to F*^^': 

(61) 



E[e-/o'-='^«|5f)=a;] ^^ 



Prom SIR-SV to DIR-LV Let us now specify a Heath Jarrow and Morton 
(1992) diffusion for the interest rate model and see precisely how one could 
extract, using Gyongy's result, the volatility of the forward contract under 
deterministic interest rates from the volatility of the forward contract under 
stochastic interest rates. Let us recall that in a standard HJM framework, the 
instantaneous forward rate follows 

d/(t, T) = f cr(t, T) I a{t, u)du\ dt + cr(t, T)dWl 

where cr(i, T) is a stochastic process adapted to its canonical filtration and where 
the price satisfies 

Bit,T)=cxp(- J^ f{t,s)di 

By definition rt — ,f{t,t) and then we obtain 
dB{t,T) 



ndt - aB{t,T)dWi 



B{t,T) 

r 

aB{t,T) = / iT{t,u)du 
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For our purpose, let us consider a general model 

^ - ndt + ./VtdWt 

then recall the price of the T-forward contract written on S 



Fl 



B{t,T) 



where we assume d < W, W^ >t= pdt. We are now able to write the dynamics 
of F^ under Q the risk-neutral measure: 

dF^ dSt d<S.,B{-,T)>t (dB{t,T) d<B{-,T)>t 



Fj St StB{t,T) \B{t,T) B^t,T) 

If we introduce the T-forward probability measure as above by 

dO^ e-.C '■='''' 



dQ B{0, T) 

we explain the dynamics of F^ under Q 

— f = ^tdWt + (TB(t, T)dWl 

where W and W^ are Brownian motions under Q-^ such as d < VF, W^ >t= pdt. 
We now apply Theorem 2.11 and obtain the existence of a process F^ solution 
of an inhomogeneous Markovian stochastic differential equation 

where /? is a Brownian motion and 

S|(i, x) = E'^[Vt + 2p^/Vt<JB{t, T) + a|(t, T)\St = xB{t, T)] 
If we consider a local volatility model with deterministic interest rates as follows 

'^=^f(0,t)dt + o{t,St)d^t 
ot 

the dynamics of the T-forward contract then becomes 
dFT 



Fl 



= a(i,F,^e-/"/(0'^)'^^)d7t 



and Gyongy's result enables us to conclude that 

Y.T{t,x) =a(i,a;e- ■''*'■/("'-)'''') 
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Hence, we have proven a new relation that Hnks a local volatility framework 
with deterministic interest rates to a stochastic volatility one with stochastic 
interest rates, namely 

a^it,x) = E'^iVt + 2pVviaBit,T) + al{t,T)\St = xB{t,T)e^f^^ f^"''^""'] (62) 

An illustration of this formula can be found for a Black and Scholes (1973) 
framework with random rates for example in Hull and White (1994). 
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6 Conclusion 

This paper recalls well-known results on local volatility and establishes links 
to stochastic volatility through the powerful theorems of Krylov and Gyongy. 
These general results are then illustrated with explicit computations of local 
volatility in different stochastic volatility models where the volatility process is 
a time-space transformation of Bessel processes. In this framework, we show 
the impact of the stock-volatility correlation on the local volatility surface. 
The local volatility extracted from a stochastic volatility model allows us to get 
a precise idea of the skew generated by a stochastic volatility model. Hence, an 
important theoretical and numerical advantage of generating a local volatility 
surface from a stochastic volatility rather than from market option prices is 
the stability and the mcaningfulncss of the surface. Indeed, the local volatility 
surface constructed with the Forward PDE equation is known to be completely 
unstable whereas as one can see the one built from a stochastic volatility is 
really smooth. 

With the growth of hybrid products, it has been important to seriously consider 
the issue of volatility calibration under stochastic interest rates and that is the 
reason why we exhibit different relations between local volatilities, stochastic 
volatilities and derivative prices. It is shown that Dupirc (1994) and Derman 
and Kani (1998) formulas can easily be extended and that it is possible to 
relate any continuous stochastic volatility model with stochastic interest rates 
to a local volatility one with deterministic interest rates. By extending the local 
volatility formula to a stochastic rates framework, it is observed that a market 
premium for the hybrid correlation risk is to be implied for the construction of 
the local volatility surface, which can be performed under the (i7C)-IIypothesis 
as at some point a market premium for the volatility risk is to be taken into 
account. 

A remaining interesting question is the existence of a local volatility diffusion 
with a general Ito interest rates process framework such that the joint law of 
instantaneous rate, the discount factor and the stock price is the same as the 
one in a stochastic volatility framework. 
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